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ABSTRACT 

High dispersion near-infrared spectra have been taken of seven highly- evolved, 
variable, intermediate-mass (4-6 M Q ) AGB stars in the LMC and SMC in order to 
look for C, N and O variations that are expected to arise from third dredge- up and 
hot-bottom burning. The pulsation of the objects has been modelled, yielding stellar 
masses, and spectral synthesis calculations have been performed in order to derive 
abundances from the observed spectra. For two stars, abundances of C, N, O, Na, Al, 
Ti, Sc and Fe were derived and compared with the abundances predicted by detailed 
AGB models. Both stars show very large N enhancements and C deficiencies. These 
results provide the first observational confirmation of the long-predicted production 
of primary nitrogen by the combination of third dredge-up and hot-bottom burning 
in intermediate-mass AGB stars. It was not possible to derive abundances for the 
remaining five stars: three were too cool to model, while another two had strong 
shocks in their atmospheres which caused strong emission to fill the line cores and 
made abundance determination impossible. The latter occurrence allows us to predict 
the pulsation phase interval during which observations should be made if successful 
abundance analysis is to be possible. 

Key words: stars: AGB and post-AGB - stars: abundances - stars: oscillations. 



1 INTRODUCTION AGB stars is a possible source for this primary nitrogen. 

Detailed stellar evolution and nucleosynthesis calculations 

AGB stars are predicted to be major contributors to the sug gest that the hot-bottom burning in AGB stars may also 

. enrichment of the interstellar medium in carbon nitre- involve the Ne _ Na and Mg _ A1 chainSj leading to changes 

gen and s-process elements (e.g. | lben and Trurarj [197&J; in the sur f ace abundances of Al am d Mg as well as ni- 

| Dray et al. | IgOOj). The 1 C and s-process elements are troge n (iKarakas and Lattanzid Efflol lLattanzio and Wood! 

brought to the surface of AGB stars during the third dredge- fioM\ . Unfortunately, all the evolution and nucleosynthe- 

up (| lben and Renzini| | l983| ) and they are then ejected into sig calculations are very uncertain due to the problem of 

the interstellar medium by the extensive stellar winds that treating convection, especially the critical process of convec- 

occur during the AGB phase of evolution. The 14 N produc- tive overshoot . Observational constraints from intermediate- 

tion is mainly from the envelopes of the more massive AGB mags AGB stars are required to constra i n the convection 

stars where the hot-bottom burning process ([Scalo et al. | parameters such as mixing length and overshoot distance, 

[1975|) can cycle the entire envelope through the outer parts which are critical to the nuc i eosynt hesis results. This study 

of the hydrogen-burning shell, converting C to N. In was initiated to provide such constraints, 
the case where the 12 C in the envelope has a component 

that comes from third dredge-up, the 14 N resulting from Considerable numbers of intermediate-mass AGB stars 

the dredged-up 12 C is of primary origin. There is strong are known in the Magellanic Clouds (Wood, Bessell & Fox 

evidence that some of the 14 N in the universe is of pri- 1983), and these stars can be used to provide the obser- 

mary origin (|van Zee et al.|[l998l ). with the primary source vational tests required for the evolution and nucleosynthe- 

dominating for metal abundances less than one third so- sis calculations. So me work has been done in this area by 

lar. Hot-bottom burning in intermediate-mass (M > 3M Q ) ISmith et al] (|2002l ) who examined the C, N and O abun- 
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dances in a sample of luminous AGB stars in the LMC. 
They found that 14 N was enhanced in a way consistent with 
first dredge-up only, without the n eed for hot-bottom burn- 
ing. However, the stars studied bv lSmith et al.l (120021 ) were 
not known variable stars, meaning that they were not near 
the end of their AGB evolution where the effects of dredge- 
up and hot-bottom burning should be most pronounced. In 
addition, the lack of pulsation means that there was no way 
to estimate the mass of the stars involved. 

In this study, we have aimed to get surface abundances 
for a small sample of intermediate-mass, pulsating AGB 
stars. These stars are near the end of their AGB evolution 
and are about to enter their final superwind stage where 
most of the material currently in their envelopes will be 
ejected back into the interstellar medium (in fact, two of 
the stars are IRAS sources already exhibiting strong super- 
winds). Since they are pulsating, we can also estimate their 
current masses from pulsation theory. On the other hand, 
because of their low effective temperatures and their pulsa- 
tion, the model atmospheres required for the interpretation 
of their spectra are complicated and difficult, and must in- 
volve the dynamics of the atmosphere. 

In Section [2] we describe the sample of stars and the 
photometry and spectra obtained for them. In Section [3] 
we describe the pulsation models, in Section [3] the model 
atmospheres are described and the derived abundances are 
given. The results are discussed in the final section. 



2 OBSERVATIONS 

2.1 Near-infrared photometry 

The sample of objects observed in this study is given in 
Table [1] All these stars are oxygen-rich, pulsating, lumi- 
nous AGB stars. All the stars have been monitored in J 
and K with the ANU 2.3m telescope at Siding Spring Ob- 
servatory, using either a single channel pho tometer (prior to 
Marc h 1994) or the infrared array CASPIR (|McGregor et all 
1 19941 ). The mean bolometric luminosity Mboi of each star 
was derived from the JK photo metry and monitor i ng, using 
the bolometric corrections from iHoudashelt" et al.1 (|2000ai r5) 
with [Fe/H] = -0.3 and -0.6 for the LMC and SMC, respec- 
tively. Distance modulii of 18.54 and 18.93 and reddenings 
E(B-V) of 0.08 and 0.12 were assumed for the LMC and 
SMC, respectively (Keller and Wood 2006). 

The least evolved of the stars in Table (TJ with lowest 
luminosity, smallest amplitude and shortest pulsation pe- 
riod, is the most luminous red giant in the LMC cluster 
NGC 1866. At the other extreme are two IRAS sources with 
very long pulsation periods and large amplitudes. The lu- 
minosities and periods of the stars in the sampl e mean that 
their current masses lie in the range 3-8 Mq (|Wood et al.l 
1 19831 ). 

2.2 The high-resolution near-infrared spectra 

Spectra were taken of these stars using the Phoenix spectro- 
graph (Hinkle et al. 2003) on the Gemini South telescope. 
The observations were done in service mode, and the obser- 
vation dates are listed in Table[T] For each star, observations 
were made in 3 separate bands, as shown in Figures [1] [2] and 



Table 1. The intermediate-mass AGB star sample 



Star 


A-/ b oi 


P 1 


AK 


MC 


Date 2 


NGC 1866 #4 


-6.00 


158 


0.10 


LMC 


030210 


HV 2576 


-6.61 


530 


0.25 


LMC 


030209 


HV 11303 


-5.59 


534 


0.70 


SMC 


030728 


HV 12149 


-6.84 


742 


0.70 


SMC 


020920 


GM 103 


-6.80 


1070 


1.30 


SMC 


030728 


IRAS 04516-6902 


-6.80 


1090 


1.30 


LMC 


030211 


IRAS 04509-6922 


-7.33 


1290 


1.50 


LMC 


021202 



1 Pulsation period in days. 

2 Date of Gemini observation in the form yymmdd. For GM 103 
and HV 11303 the observations were spread over 6 days and the 
given date is the mean value. 



Table 2. The modelled stars 



Star 


L/L Q 


T e ff 


M/M Q 


e/u p 




NGC 1866 #4 


20000 


3490 


4.0 


2.57 


0.05 


HV 2576 


35000 


3350 


6.0 


2.21 


0.51 


HV 11303 


34000 


3490 


4.6 


2.21 


0.58 


GM 103 


41840 


3040 


6.0 


1.40 


1.30 



Turbulent viscosity parameter. 



[3] Firstly, an observation centred near 1.554 was made 
so that the OH lines could be measured for the derivation 
of an oxygen abundance; lines of CN also occur in this piece 
of spectrum. Secondly, an observation near 2.340 pm was 
made to measure 12 CO (and possible 13 CO) lines and hence 
the 12 C abundance; given the C abundance, the CN lines 
near 1.554 fim then yield a 14 N abundance. The Na abun- 
dance was also obtained from this spectral region. Thirdly, 
an observation at 2.112 fj,m was made to include lines of Al 
(and possible Mg). Lines of Fe, Sc and Ti also occur in the 
spectra, giving an estimate of the metal abundance. 



3 PULSATION MODELS 

Pulsation models and model atmospheres were made for 
four stars in Tabled NGC 1866 #4, HV2576, HV 11303 
and GM 103. The spectra of the remaining stars in Table 
[T] (HV 12149, IRAS 04516-6902 and IRAS 04509-6922) were 
considered too difficult to model at this stage because of 
the extreme appearance of the spectra, with no clear con- 
tinuum remaining (see Figs. [T][2]and[3]|. This is presumably 
a result of the very large amplitude and very low effective 
temperatures of these stars. The J and K photometry for 
the four stars for which modelling was attempted is given in 
Table [3] Typical photometric errors are less than 0.03 mag- 
nitudes. The J and K magn itudes have been conver ted to 
those on the AA O system of I Allen and Cragd |l98St ) using 
the conversions in iMcGregorl (|l994l ). 

For each observed star, given estimates of the mean lu- 
minosity L and the effective temperature T ff from the J and 
K photometry, the stellar radius was computed from the 
definition L = 4naR 2 T^ H . Then, given the radius and the 
known pulsation period, the current stellar mass was com- 
puted from the P-M-R relation of stellar pulsation. These 
values of mass and luminosity were then used in the nonlin- 
ear pulsation calculations for the star. 
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Figure 1. The spectra of all objects in the region near 1.554 fim. The most important lines used for abundance determination 
marked. 
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Figure 2. The same as Figurc[T]but for the region near 2.112 fim. 
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Figure 3. The same as Figurc[T]but for the region near 2.340 /an. 



Table 3. JK photometry 



JD24 1 


J-K 


K 


JD24 


J-K 


K 


NGC1866 #4 




HV 11303 (continued) 


50731 


1.23 


9.71 


52804 


1 


29 


10.02 


50772 


1.27 


9.69 


52924 


1 


26 


9.97 


50885 


1.24 


9.68 


53516 


1 


20 


9.71 


51157 


1.30 


9.67 


53574 


1 


20 


9.38 


52564 


1.29 


9.71 










52714 


1.33 


9.64 


GM 103 








52803 


1.25 


9.69 


45980 


1 


66 


9.55 


52924 


1.28 


9.69 


46341 


1 


19 


8.65 








46645 


1 


54 


8.57 


HV2576 2 






48142 


1 


63 


9.31 


52563 


1.39 


9.00 


48168 


1 


70 


9.23 


52715 


1.43 


9.13 


48640 


1 


64 


8.36 


52804 


1.29 


8.97 


48851 


1 


62 


8.78 


52924 


1.36 


8.87 


48934 


1 


63 


9.13 








48992 


1 


67 


9.33 


HV11303 2 - 3 




49259 


1 


89 


9.80 


49259 


1.19 


9.81 


49317 


1 


87 


9.82 


49291 




9.61 


49374 


1 


68 


9.45 


49317 


1.17 


9.47 


52564 


1 


65 


9.33 


49374 


1.19 


9.32 


52714 


1 


59 


8.74 


52564 


1.21 


9.27 


52803 


1 


61 


8.61 


52715 


1.39 


9.56 


52924 


1 


64 


8.40 



1 JD24 is Julian Date - 2400000. 2 Extra photometry is given in 
Wood ct al. (1983). 3 Extra photometry is given in Catchpole and 



Feast (1981). 



The nonlinear pulsation code described in Keller and 
Wood (2006), and references therein, was used in this study. 
The turbulent viscosity parameter a v was adjusted to give 
the observed K light curve amplitude. The K and V mag- 
nitudes were computed using L and T e g from the pulsa- 
tion models together wit h the bolometric corrections from 
iHoudashelt et all (j2000al lbh. 

For three of the four modelled stars, MACHO light 
curves are available: the V magnitude was computed from 
the MACHO magnitude using the transforms in Bessell and 
Germany (1999). For the fourth star, HV 11303, an Eros2 
light curve was kindly provided by Patrick Tisserand (pri- 
vate communication), who also provided a conversion from 
the two Eros2 bands to the V band. The combined V and 
K light curves give a V-K colour curve. The effective tem- 
peratures of the models were adjusted (by altering the ra- 
tio of mixing-length to pressure scale height ^/H p to repro- 
duce the observed V-K colour when available, or the J-K in 
the one remaining case of HV 11303 (the Eros2 light curve 
was obtained after the modelling was completed). The final 
model parameters are given in Table [2] In all CEIS6S, £L helium 
abundance Y = 0.30 was assumed, while a metal abundance 
Z=0.01 was assumed for the LMC stars and Z=0.004 was 
assumed for the SMC stars. 

In order to construct a model atmosphere taking 
account of the dynamical processes in th e outer lay- 
ers, p r ocedures similar to th o se described in Bessell et al.l 
Jl99rJ).[Hofmann et all ljl99ct ). IScholz and Wood! (|2000l ) and 
I Ireland et al. ( 2004al lbTl were used, with u pdate s to the 
molecular opacity from llreland and Scholzl ((20061 ). Briefly, 
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a pulsation model at the same pulsation phase as the phase 
of the Gemini observation was extracted from each pulsa- 
tion series. The radial profiles of gas pressure and velocity 
for this model were transferred to a non-grey, spherical at- 
mosphere code, and the temperature structure was iterated 
to a new state, now consistent with the atmosphere code. 
Once the temperature, gas pressure, density and velocity 
structure were all consistent in the atmosphere code, the 
model atmosphere structure was transferred to a separate 
spectrum synthesis code which was used to compute line 
profiles (see Section [!}. 

In the atmos phere code, half - and quarter-solar metal 
abundances from iGrevesse et al. I (|l996h were adopted for 
LMC and SMC stars, respectively. Radiative and local ther- 
modynamic equilibrium was adopted and, in particular, the 
shock-heated region behind the shock front was assumed to 
have negligible width and to have no effects on the temper- 
ature stratification. 

The pulsation series and the extracted dynamical model 
used as input to the atmosphere code is now described for 
each star. 



3.1 NGC 1866 #4 

NGC 1866 #4 is the most luminous red giant in the LMC 
cluster NGC 1866 (the designation #4 is that given by Fro- 
gel, Mould and Blanco 1990: note that Maceroni et al. 2002 
have mislabelled #4 as #2 in their paper). The MACHO 
light curve of this star shows a period of 158 days . Mod- 
elling the cluster HR-diagram, iBrocato et all (|2003T l derive 
masses for the He-burning giants of 3.8-4.15 Mq. Our de- 
rived mass for NGC 1866 #4 is 4 M . With this mass and 
the computed mean luminosity and T e g, this star must be 
pulsating in the first overtone mode, a result also consistent 
with the small amplitude. 

Fig.f3]shows the observed and computed K, V and M^oi 
light curves. The period and amplitude of the pulsation and 
the V-K colour are reproduced well. The Gemini observa- 
tion was made at a phase of ~0.1 after visual maximum. 
The structure of the pulsation model at this time is shown 
in Fig. [5] The pulsation amplitude of this star is quite small 
so there are no shock waves present in the atmosphere. The 
temperature solution from the model atmosphere program 
produces a cooling of the outer layers and a warming at 
moderate optical depths compared to the modified Edding- 
ton approximation used in the pulsation models. 



-6.2 
| -6.1 

a 

-6 

-5.9 

14.5 

> 15 

15.5 

16 
9.5 

9.6 

9.7 



NGC 1866 #4 M=4M Q L^OOOOLg Z = 0.01 



xxxx„x x 



observation 



-•— I- 



0.5 1 
Phase 



1.5 



Figure 4. The K, V and M bol light curves of NGC 1866 #4. 
The solid points are model values (several cycles are overplotted) , 
while the crosses are observations or M\> \ values computed from 
the J and K photometry in Table [3] The V values are computed 
from the MACHO Mg and Mr photometry using the transforms 
of Bessell and Germany (1999). The arrows show the phase of 
the Gemini observations and the phase of the pulsation model 
extracted for model atmosphere computation. 



period and amplitude of the pulsation and the large V-K 
colour (V-K ~ 7) are reproduced well. The Gemini observa- 
tion was made near minimum light of the K light curve, or 
between the double-humped maximum of the optical light 
curve. The structure of the pulsation model at this time 
is shown in Fig. [7] The velocity gradient through the atmo- 
sphere at the time of observation was small and there was no 
shock wave present in the atmosphere. As with NGC 1866 
#4, the temperature solution from the model atmosphere 
program produces a cooling of the outer layers and a warm- 
ing at moderate optical depths compared to the modified 
Eddington approximation used in the pulsation models. 



3.2 HV2576 

HV2576 is an LMC star which is more luminous than 
NGC 1866 #4. It has a much longer period, a larger ampli- 
tude, and it pulsates in the fundamental mode. The derived 
mass for HV2576 is 6 M Q . 

Fig. [6]shows the observed and computed K, V and Afboi 
light curves. Like many of the intermediate-mass, large- 
amplitude AGB pulsators in the Magellanic Clouds, this star 
shows a double-humped optical light curve. The pulsation 
models do not accurately reproduce the double hump but 
they do show some evidence for it. The K light curve does 
not show evidence for a prominent hump, but the points are 
too sparse in the light curve to be definitive. Overall, the 



3.3 HV 11303 

HV 11303 is an SMC star which is almost a twin of HV 2576 
in the LMC in terms of its pulsation period and luminosity. 
However, it has a much larger pulsation amplitude. It pul- 
sates in the fundamental mode and the derived mass is 4.6 
Mq. This mass is smaller than that of HV2576, which may 
explain the larger pulsation amplitude. 

Fig. |S]shows the observed and computed K, V and Mboi 
light curves. The pulsation models reproduce the observed 
K light curve reasonably well. The Gemini observation was 
made near minimum light of the K light curve. The structure 
of the pulsation model at this time is shown in Fig. [9] Note 
that there is an uncertainty of about 0.1 in the phasing of 
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380 400 420 



440 



Figure 5. The density p, temperature T and velocity v plotted 
against radius R in a pulsation model of NGC 1866 #4 at the 
same phase as the phase of Gemini observation of the star. The 
dashed temperature profile is that obtained from the non-grey 
model atmosphere code. 
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Figure 7. The same as Fig.|5]but for HV2576. 
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Figure 8. The same as Fig.g]but for HV 11303. 



Figure 6. The same as Fig. E| but for HV2576. 
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Figure 9. The same as Fig. E] but for HV 11303. 



Figure 10. The same as Fig.Q]but for GM103. 



the models relative to the K light curve. With the adopted 
phasing, a strong shock was just emerging into the atmo- 
sphere at the time of Gemini observation. In addition, there 
was a large velocity gradient through the atmosphere above 
the shock. The temperature solution from the model atmo- 
sphere program shows a large temperature drop in the outer 
layers. This behaviour is typical for extended atmospheres 
in a certain parameter range where stron g water abso rption 
abruptly dominates in the high layers (|Scholdll985h . The 
warmer, grey pulsation model used as input to the atmo- 
sphere code does not have water-forming low temperatures. 



3.4 GM103 

GM 103 (Groenewegen et al. 1995) is also an SMC star but 
of greater luminosity, amplitude and period than HV 11303. 
It pulsates in the fundamental mode and the derived mass 
is 6 M . 

Fig. [TO] shows the observed and computed K, V and 
Aiboi light curves. The pulsation models reproduce the K 
light curve and the large V-K colour (V-K ~ 8) reasonably 
well. However, the V light curve of the pulsation models 
seems to rise towards maximum earlier than the observed 
light curve. The Gemini observation was made near maxi- 
mum visible light, or about phase 0.1 before maximum of 
the K light curve. The structure of the pulsation model at 
this time is shown in Fig. 1111 This star has a strong shock 
situated in the middle part of the atmosphere at the time of 
Gemini observation, with moderate velocity gradients above 
and below the shock. The temperature solution for GM 103 
from the model atmosphere program shows a larger devi- 
ation from the pulsation solution than in any other star. 
This is due to the strong density drop at the shock, with 
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Figure 11. The same as Fig. El but for GM103. 



© 2007 RAS, MNRAS OOO.IlMil 



8 J. A. McSaveney et al. 



the model atmosphere's water-dominated temperature drop- 
ping rapidly near what is the effective edge of the star at the 
shock front. 



4 SPECTRAL SYNTHESIS AND 
ABUNDANCE DERIVATION 

4.1 Line synthesis calculations 

The spherical model atmospheres described in Section [3] 
were put into a compatible spectrum synthesis program (us- 
ing spherical geometry) to generate line spe ctra. The spec - 
trum synthe si s pro gram was that used by ISchobJ (| 19921 ). 
iBessell et ail (|l996ft and IScholz and Wood! (|2000l ). except 
that partial pressures of a larger number of molecules were 
calculated using a computer program supplied by K. Ohnaka 
(private communication). This progr am use s an e quation of 
state based on molecular constants oflTsuiil dl973l l , updated 
for CN using da ta from ICostes et all (|l99Ct ) and for TiO 
from iTsuiil (|l978l ). Lines are assumed to be formed in local 
thermodynamic equilibrium. As shock-heating behind the 



shock front is neglected in these models, emission compo- 
nents of synthetic line profiles may only occur as a conse- 
quence of large atmospheric extension (P-Cygni-like emis- 
sion; see Scholz 1992). 

The overall spectrum in each of the three filter re- 
gions was computed using a reasonably comprehensive set 
of atomic and diatomic molecular lines, although only 
specific lines were used in the actual abundance deriva- 
tions. For atomic lines, wavelengths, excitation potentials 
and log gf values were taken from the V ienna Atomic 
Line Database (VALD, iKupka et all Il999h , with param- 
eters of lines selected for abundance determ ination be- 
ing checked against those of lSmith et alj (|2002l ). Molecular 
line data proved more problem atic as no s i ngle s ource ex- 
ists. CO lines were taken from jGoorvitchl (|l994h and CN 
lines from Aoki and TsiiiH (1 1997 ) and the SCAN database 
(|jorgensen and Larssonl 199Cj[). CH lines were also obtained 
from the SCAN database ( Jorgensen et al.l Il996h . In each 



of these cases, line s were cross-checked against those from 
ISmith et all | |2002ft . The OH lines proved most problem- 
atic. _Aji_Jncomplete sel ection of lines was ass embled from 



ISmith et all (|2002h and iMelendez et all d200ll) : seve ral ad- 
ditional lines listed in the source Abrams et al.l (|l994h could 
not be not included due to a lack of gf values. 



4.2 The derivation of abundances 

Starting abundances were adopted by approximating LMC 
and SMC metal abundances as half and quarter sola r, re- 
spectively, with solar defined from lAsplund et alj (|2005l ) (see 
Table [5}. The CO and OH lines computed with these abun- 
dances were examined to establish the microturbulent ve- 
locity. 

As these relatively massive stars are expected to have 
altered CNO abundances as a result of nucleosynthetic pro- 
cessing, a serie s of CNO values base d on the AGB evolution 
calcul ations of iKarakaj (|2003T ) and iKarakas and Lattanziol 
(|2003h was adopted, and the resultant spectra were then 
calculated. The CNO combination giving the best overall 
fit to the observed spectra was used as the primary starting 
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Figure 12. Fitted OH lines for HV2576. In this and subsequent 
figures, fits are shown for the three abundance values shown on 
the figure (where A(X) = log[n(X)/n(-ff)] + 12). The dashed lines 
correspond to changes of ±0.15 dex from the adopted abundance 
(solid line). 
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Figure 13. Fitted OH lines for NGC 1866 #4. 

point for subsequent refinement of the individual CNO abun- 
dances. Once the best CNO abundances had been found, the 
metallic abundances were refined to give a best fit. Lines 
used for final abundance derivations are listed in Table [4] 
and a selection of the fitted lines are shown in Figs. [12] to 

m 

The abundances derived from different lines of the same 
species show some scatter, with the adopted lines generally 
found to yield reasonably consistent values. Figs. 1141 and 1151 
show examples of the fit errors for the CO lines. It is clear 
that neither CO line is fitted exactly by the final adopted 
abundance (that of the mean of the best fits to 6 lines), but 
both lines are fitted reasonably well, with one too strong and 
one too weak. Figs. PT21 and 1131 where Fe I lines are present 
on either side of the OH lines, show further examples of 
fit errors. The adopted Fe abundance is based on a fit to 
the 2.1095^m line in Figs. [24] and [25] However, in HV2576 
(Fig. [12]), this abundance is too high for the 1.553424/xm Fe 
line to the blue and too low for the 1.553757/im line to the 
red (or, more likely, the models are missing an unidentified 
line at this wavelength). For NGC 1866 #4 (Fig. [13]), the 
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Figure 14. Fitted CO lines for HV2576. 
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Figure 15. Fitted CO lines for NGC 1866 #4. 

1.553424/xm line fits well, while once again the 1.553757/im 
line is too weak (or an unidentified line is missing in the 
model) . 

There were also cases where observed and computed 
line velocities do not match to better than a few km s _1 , as 
can be seen in Figs. 1141 to [TT1 The velocity mismatch could 
be due to a number of potential causes, including inexact 
adopted line wavelengths, inexact wavelength calibration, or 
the model velocities as a function of depth not being correct. 

The microturbulent velocity required for line broad- 
ening varied from star to star. The adopted value for 
NGC 1866 #4 (3 km s" 1 ) is close to the value used by 
ISmith et al.l (|2002h . However, a larger value of 7 km s _1 was 
required for HV2576. This value was derived by fitting to 
the CO lines. Slightly different values would be obtained by 
fitting different lines. This is probably a reflection of different 
amounts of turbulence at different levels in the atmosphere: 
the CO lines tend to be broader than the metallic lines and 
they are, on average, formed further out in cooler layers. 

Our final derived abundances are shown in Table [5] For 
the elements Fe, Na, Al, Sc, Ti, and N, the error estimates 
in the table are those arising from the fit to the line only, as 
only one line was available. They are eye-estimates based on 
the fits 0.15 dex above and below the adopted optimum fit 
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Figure 16. Fitted CN line for HV2576. 
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Figure 17. Fitted CN line for NGC 1866 #4. 
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Figure 18. Fitted Sc I line for HV2576. 
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Figure 19. Fitted Sc I lino for NGC 1866 #4. 



Figure 22. Fitted Na I line for HV2576. 
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Figure 20. Fitted Ti I line for HV2576. 



Figure 23. Fitted Na I line for NGC 1866 #4. 
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Figure 21. Fitted Ti I line for NGC 1866 #4. 



Figure 24. Fitted Fe I and Al lines for HV2576. 
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Figure 25. Fitted Fe I and Al lines for NGC 1866 #4. 

Table 4. Spectral line data. 
A x Source 

VALD 

Smith et al. (2002) 



(A) 




(eV) 


i°g(s/) 


Fe I 








21095 


446 


6.20 


-0.69 


Na I 








23378 


945 


3.75 


-0.420 


Sc I 








23404 


756 


1.44 


-1.278 


Ti I 








15543 


720 


1.88 


-1.481 


Al I 








21093 


029 


4.09 


-0.31 


OH 








15535 


489 


0.51 


-5.23 


15536 


707 


0.51 


-5.23 


15560 


271 


0.30 


-5.31 


15565 


815 


0.90 


-5.00 


CO 








23396 


260 


0.37 


-5.17 


23398 


224 


1.73 


-4.44 


23406 


345 


0.00 


-6.57 


23408 


509 


0.36 


-5.19 


23424 


328 


1.80 


-4.427 


23426 


322 


0.00 


-6.565 


CN 








15563 


367 


1.15 


-1.141 



Smith et al. (2002) 
Smith et al. ('20021 



VALD 



Smith ct ah (20021 
Melendez et al. (2001) 
Smith et al. (20021 
Melendez et al. (2001) 



Smitli^et^L 
Smith ct al. 
Smith_et i ^l 1 
Smith ct al. 
^mitli^et^L 
Smith et al. 



[2002) 
f2002) 
f2002) 
f2002) 
[2002) 
f2002) 



Smith et al. (2002 



abundance. The C and O values are based on averages of the 
fits to 6 and 4 lines, respectively, with the stated errors being 
based on the standard deviation. The lines used are given in 
Table [4] Systematic errors, such as errors in the pulsation 
model and its phase, will undoubtedly add significantly to 
the fit error, but such additional errors are hard to estimate. 



4.3 HV 11303 and GM 103 

As noted in Section^ the SMC stars HV 11303 and GM 103 
both had strong shock fronts in their atmospheres at the 
time of observation with Gemini South. This has made 
abundance derivation with our atmosphere code impossible. 



Fig. [26] shows a sample of the spectrum for each star, along 
with an attempt at synthesisizing the spectrum using the 
pulsation models from Section [3] 

In the case of GM 103, the synthesized absorption lines 
are much stronger than the observed lines. A comparison of 
the observed and synthesisized spectra shows clearly that 
the observed line cores are filled in by strong emission in ev- 
ery case. An alternative explanation could be that each line 
is split into a pre- and post-shock absorption component, 
with no emission in the middle. However, the weakness of 
the observed absorption lines compared to the model lines 
suggests that there is indeed a strong emission component. 
In HV 11303, the lines are much broader than in GM 103, 
with evidence for an emission component as well as absorp- 
tion components corresponding to the pre- and post-shock 
velocities. 

As noted in Section |3j the spectral synthesis code ig- 
nores the immediate post-shock, high-temperature region 
so that there is no post-shock emission component in our 
synthesized spectra. When a strong shock is present in the 
atmosphere, this model deficiency prevents abundance de- 
termination from lines which have a post-shock emission 
contribution. 

These results show that observations useful for abun- 
dance analysis must be obtained at phases of the pulsation 
cycle when no strong shock wave is present in the atmo- 
sphere. In these relatively massive stars, the shock wave en- 
ters the lower part of the stellar atmosphere near minimum 
of the K light curve, or phase ~0.2-0.3 past minimum vi- 
sual light. It has passed through the line forming part of 
the atmosphere half a pulsation cycle later, corresponding 
to ~0.1-0.2 in phase past K maximum or after phase ~0.4 
of the visual light curve. 



5 DISCUSSION 

The abundances obtained here for NGC 1866 #4 and 
HV2576 are shown with other comparable observations in 
Fig. 1271 The Fe abundance i s close to hal f solar , consistent 
with t he values obtained by ISmith et all (|2002l ). iHm" et al] 
(|200d ) and other studies of young objects in the LMC. The 
other metal abundances are similarly near half-solar, except 
for Ti in HV2576 where our derived value seems anoma- 
lously low. 

The CNO abundances are the most interesting part of 
this study, since in a highly-evolved, intermediate-mass AGB 
star they will be altered by the 1st, 2nd and 3rd dredge-ups 
and by hot-bottom burning. 

Evidence for hot-bottom burning has been found in 
intermediate-ma ss A GB stars in the Ma g ellanic Clouds by 
iPlez et al] (| 19931 ) and lSmith et~ai1 (|l995l ). |Plez et al] (|l993T ) 
examined spectra of 7 luminous AGB stars in the SMC and 
found them to be Li-rich, C-poor and with l ow 12 C/ 13 C 
ratios , consistent with hot-bottom burning. ISmith et al] 
\ 19951 ) examined a large number of SMC and LMC stars 
and found many luminous AGB stars, including HV2576, 
showed strong Li lines that could be attributed to hot- 
bottom burning. In fact, HV2576 had the l argest estimated 
Li ab undance of all the stars examined. iMaceroni et al] 
(|2002l ) examined Li lines in the giant stars in NGC 1866 
and found that NGC 1866 #4 (which they mislabel as #2) 
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Figure 26. Top: A section of the 2.34 /im spectrum of GM 103 (points), overlayed with a synthesised spectrum (continuous line). Bottom: 
The same piece of spectrum for HV 11303. 



Table 5. Derived abundances 



Star 



A(Fe) 



A( 12 C) 



A( 14 N) 



A( 16 Q) 



A(Na) 



A(Sc) 



A(Ti) 



A(A1) 



LMC 
HV 2576 
NGC 1866 #4 
Smith et al. 2002 

Half-solar 



7.20 ± 0.10 7.01 ± 0.20 8.50 ± 0.05 8.28 ± 0.08 5.35 ± 0.10 2.90 ± 0.10 3.60 ± 0.05 5.50 ± 0.05 
7.05 ± 0.15 7.36 ± 0.16 8.55 ± 0.15 8.80 ± 0.17 5.75 ± 0.15 3.40 ± 0.15 4.70 ± 0.15 5.75 ± 0.10 
6.37 - 7.16 6.53 - 7.86 7.14 - 8.24 7.82 - 8.33 4.69 - 5.84 2.01 - 2.91 3.84 - 4.61 



7.15 



8.09 



7.48 



,36 



5.87 



2.75 



4.60 



5.87 



Notes: A(X) = log[n(X) /n(H)] + 12. Half-solar abundances are half the solar abundances in lAsplund et aD ||2005T ). The errors are 
estimates of the line fitting error alone - see Sect. 11.21 



had enhanced Li compared to other stars in the cluster, in- 
dicating that hot-bottom burning is operating in this star 
also. 

There is evidence that the third dredge-up has been 
operating i n HV 2 576, as well as hot-bottom burning. 
IWood et afl (|l98ot ) give a spectral type M5.5/S1 for 
HV2576, indicating the presence of excess ZrO in the spec- 
trum, which is a result of t he dredge-up of s- process ele- 
ments during thermal pulses. ISmith et al.l (|l995l ) noted that 
HV 2576 had one of the strongest ZrO bands in their large 
sample of stars, indicating that a large amount of third 
dredge-up has occurred in this star. 

The combination of third dredge-up and hot-bottom 
burning is predicted to be a significant source of primary 
nitrogen in the universe. The 12 C dredged up at a ther- 
mal pulse is converted into 14 N during the next interpulse 
period when hot-bottom burning operates during the quies- 



cent H-s h ell bu rning phase. lEattanzio and Wood! (|2004t ) and 
iKarakas] (|2003l ) (see their Figs. 2.37 and 5.6, respectively) 
show the effects of third dredge-up and hot-bottom burn- 
ing in intermediate-mass, luminous and sub-solar metallic- 
ity AGB stars such as those we are dealing with here. Ini- 
tially on the AGB, just at the onset of thermal pulses, hot- 
bottom burning by the CN cycle reduces the 12 C abundance 
to ~1/15 of its initial main-sequence value, and increases 
the 14 N abundance to ~5-6 times its initial value. The ratio 
12 C/ 14 N at this time is ~1/15. During subsequent evolution, 
the dredge-up of 12 C by third dredge- up at thermal pulses, 
followed by hot-bottom burning, causes a steady rise over 
many pulsation cycles in both the C and N abundance, with 
the C/N ratio remaining almost constant at ~1/15. Dur- 
ing this relatively long evolutionary process, there is a small 
decrease in the O abundance, but this effect would be too 
small for us to reliably detect. 
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Looking at Fig, 1271 and Table[5j we see that the C abun- 
dance in both NGC 1866 #4 and HV2576 is about 0.1 times 
half-solar. At the same time, the N abundance is ~10 times 
the half-solar value. Both these numbers are consistent with 
the strong operation of hot-bottom burning in the envelopes 
of these stars. Furthermore, the very large N abundances can 
only be explained by hot-bottom processing of dredged-up 
12 C into 14 N i.e. a large fraction of the nitrogen in the en- 
velopes of these stars must be of primary origin. This is the 
first direct demonstration that this long-predicted source of 
primary nitrogen does exist. 

The O abundance in HV2576 is essentially unchanged 
from the half-solar value, as predicted by models. There 
is some evidence for an increase in the O abundance in 
NGC 1866 #4. If confirmed, this would lend support to con- 
vective theories such as that of iHerwig et af] (|l997h which 
lead to overshoot at convective boundaries, enrichment of 
the intershell region in both 12 C and 16 at each helium 
shell flash, and subsequent dredge-up of both 12 C and le O 
in third dredge-up events. 

A strong disagreement between the observations and 
evolutionary models without convective overshoot relates 
to the total stellar mass required for the efficient opera- 
tion of hot-bottom burning. The mass of NGC 1866 #4 must 
be no more than ~4 Mq because of cluster membership. 
This star shows evidence for efficient hot-bottom burning, 
while models without convective overshoot do not predict ef- 
ficient hot-bottom burning at masses this low. This suggests 
that overshoot (inwards) of convective material at convective 
boundaries is required to significantly extends the convec- 
tive envelope i nto n uclear-burning regions. The models of 
IVentura et all (|2002l ) use a convective theory that does pro- 
duce overshoot, and their 4 Mq, half-solar metal abundance 
AGB models do have the efficient hot-bottom burning ob- 
served in NGC 1866 #4. At 6 M , both the observations of 
HV2576 and models with or without convective overshoot 
show efficient hot-bottom burning. 



The elements Li 
evolutionary models 



Na and Al c a n also be compared with 
ith et all dl995l) give A(Li) = 3.8 
for HV 2576 while iMaceroni et alj (I2002T ) give A (Li) = 1.5 
(note that these authors used static model atmospheres, and 
observations at random phases where shocks may have been 
present in the stellar atmosphere) . Models of AGB stars with 
efficient hot-bottom burning generally show an initial large 
increase in the surface Li abundance, after which it slowly 
decreases to low values (A(Li) < 1.0) while the CN cycle 
conve rts C into N - see iKarakasl ([2003) and IVentura et al.l 
(|2002l ). The high Li value derived for HV 2576, together with 
the observed high 14 N and low 12 C abundances, is consistent 
with model predictions for a 6 Mq star about one quarter of 
the way through its TPAGB phase a nd undergoing e fficient 
hot-bottom burning - see Fig. C27 of lKarakasI ( |2003D . 

The Na and Al abundances found in this study are 
marginally below the expected half-solar values and within 
the e xpected range of observed variation. Models of AGB 
stars (Karakas 2003) with efficient hot-bottom burning show 
a small increase in Al abundance and small changes in the 
total Mg abundance but these are too small to detect here. 
The 25 Mg/ 24 Mg ratio can increase greatly in such stars, but 
we are unable to detect isotopic changes. 



a Solar (Asplund et al 2005) 
o Solar -0.3 (LMC 1st approximation) 
A HV2576 
+ NGC1866 4 
* NGC1866 (Hill et al 2000) 
Red Giant range (Smith et al 2002) 



,N 



,Fe ,Na Al 



,Sc 



01 23456789 
Figure 27. Final abundances. 

6 SUMMARY 

High-dispersion, near-infrared spectra have been used to de- 
rive abundances of luminous, intermediate-mass AGB stars 
in the Magellanic Clouds. The C and N abundances derived 
here provide the first observational demonstration that third 
dredge-up at thermal pulses followed by hot-bottom burn- 
ing can produce significant amounts of primary nitrogen 
in intermediate-mass AGB stars. The observed spectra and 
modelling show that only spectra taken during the part of 
the pulsation cycle near minimum visual light, equivalent to 
the latter part of the decline of the K light curve, are useful 
for abundance analysis. 

The results in this paper come from exploratory obser- 
vations and calculations and, in the end, only two stars were 
suitable for abundance analysis. It is clear that follow-up 
work with a larger sample of objects is required to confirm 
the results found here, and to seek C and N abundances as 
a function of stellar mass, luminosity and metal abundance. 
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